Infinitely-fast diffusion in Single-File Systems 
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We have used Dynamic Monte Carlo(DMC) methods and analytical techniques to analyze Single- 
File Systems for which diffusion is infinitely-fast. We have simplified the Master Equation removing 
the fast reactions and we have introduced a DMC algorithm for infinitely-fast diffusion. The DMC 
method for fast diffusion give similar results as the standard DMC with high diffusion rates. We have 
investigated the influence of characteristic parameters, such as pipe length, adsorption, desorption 
and conversion rate constants on the steady-state properties of Single-File Systems with a reaction, 
looking at cases when all the sites are reactive and when only some of them are reactive. We find 
that the effect of fast diffusion on single-file properties of the system is absent even when diffusion 
is infinitely-fast. Diffusion is not important in these systems. Smaller systems are less reactive and 
the occupancy profiles for infinitely-long systems show an exponential behavior. 

PACS numbers: 02.70.Uu, 02.60.-x, 05.50.+q, 07.05.Tp 
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INTRODUCTION 



In onc-dimcnsional systems such as zeolites or other 
porous structures, diffusion is a very important process. 
The pores of these structures that have the cross section 
somewhat larger than a benzene molecule, are modelled 
by Single-File Systems. In these systems particles move 
in a concerted fashion, as they are unable to cross each 
other. This process of Single-File diffusion has different 
characteristics from ordinary diffusion which affects the 
nature of both transport and conversion by chemical re- 
actions. In [0 and Q we have studied the steady-state 
and transient properties of this system. We have anal- 
ysed different situations for diffusion rates and we have 
compared the results obtained from simulation and ana- 
lytical techniques. Often diffusion is a very fast process 
compared to the other reactions in the system. We are 
thus interested to be able to model correctly infinitely- 
fast diffusion. For this purpose, we used DMC methods 
with high regular diffusion rates, assuming that these 
rates are high enough to model infinitely-fast diffusion. 

Dynamic Monte Carlo methods for very high rates are 
not very efficient and the progress of the simulation is 
slow. Moreover, considering regular reactions rates it is 
always a problem to balance between making the diffu- 
sion rates high enough so that the infinitely-fast diffusion 
effects are correctly modelled and the performance of the 
simulation. We derive here a new method to simulate 
infinitely fast-diffusion in Single-File Systems, starting 
from the Master Equation. 

The rate equations of some special limiting cases and 
an analytical description for the productivity of the sys- 
tem are also derived. We study also how the system be- 



havior changes for different sets of kinetic parameters and 
different distributions of the reactive sites. We categorize 
also interesting results obtained for profile occupancies 
for different reactive system and different distribution of 
the reactive sites. 

In section [n] we specify our mathematical model with 
the theoretical background for the analytical and simu- 
lation results. We i ntrod uce the Master Equation of the 
systems in section II B and then we simplify the Mas- 
ter Equati on rem oving the fa st reactions in section II C 
and 



In section II D 1 



II D 2 we present the simulation 



methods and we present a DMC algorithm for infinitely- 
fast diffusion. Different analytical results are presented 
In section 



III 



IV A and IV B we analyze different 



simulation results for the case when all the sites are re- 
active and when only some of the sites are reactive. We 
pay special attention to the influence of the length of the 
pipe and reaction rate constant on the site occupancy of 
the system. 



II. 



THEORY 



In this section we will give the theoretical background 
for our analytical and simulation results. First we will 
specify our model and we will derive a finite set of ex- 
act rate equations starting from the Master Equation || . 
These rate equations are used in order to derive expres- 
sions for the productivity in the system for special cases. 
We show that we can simplify the Master Equation de- 
scribing the evolution of the system over time removing 
fast reactions. We use Dynamic Monte Carlo method 
for our simulation results and we give the description of 
a Dynamic Monte Carlo-like algorithm for infinitely-fast 
diffusion. 
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A. The Model 

We model a Single-File System by a one-dimensional 
array of sites, each possibly occupied by an adsorbate. 



This is the model of diffusion and reaction in a one- 
dimensional arrangements of particles with hard-core in- 
teraction. The sites are numbered 1,2, ... ,S. A particle 
can only move to the right or to the left if an adjacent 
site is vacant. The sites could be reactive and unreactive 
and we note with JVp rot the number of reactive sites. A 
reactive site is the only place where a conversion may 
take place. 

We consider two types of adsorbates, A and B, in 
our model and we denote with Y the site occupation 
of a site, Y=(*, A, B), which stands for an vacant 
site, a site occupied by A, or a site occupied by a B, 
respectively. We restrict ourselves to the following mono 
and bi-molecular transitions. 

a) Adsorption and desorption 

Adsorption and desorption take place only at the 
two marginal sites i.e., the left and rightmost sites at 
the ends of the system. 

A(gas) + * m — ► A m 
A m — > A(gas) + * m 
B m — > B(gas) + * m , 

where subscripts m denotes a marginal site. Note that 
there is no B adsorption. B's can only be formed by a 
conversion. 

b) Diffusion 

In the pipe, particles are allowed to diffuse via hopping 
to vacant nearest neighbor sites. 

A„ + * n +l * ► *n + A n+ i 

B„ + * n +l * ► *n + B„ + j 

where the subscripts are site indices: n—1, 2, . . . , S-l. 

c) Conversion 

An A can transform into a B at a reactive site. 



By a, (3, we will indicate a particular configuration of 
the system i.e., a particular way to distribute adsorbates 
over all the sites. P a {t) will indicate the probability of 
finding the system in configuration a at time t and W a p is 
the rate constant of the reaction changing configuration 
(3 to configuration a. 

The probability of the system being in configuration 
a at time t + dt can be expressed as the sum of two 
terms. The first term is the probability to find the sys- 
tem already in configuration a at time t multiplied by 
the probability to stay in this configuration during dt. 
The second term is the probability to find the system in 
some other configuration (3 at time t multiplied by the 
probability to go from (3 to a during dt. 

P a (t+dt) = {l-dtJ2W 0a )P a (t)+dtJ2W a pPp(t) (1) 

P 

By taking the limit dt — > this equation reduces to a 
Master Equation: 



-# i =E[W a /jJ , /j(t)-W /9a P a (i)]. (2) 

at p 

Analytical results can be derived as follow. The value 
of a property X is a weighted average over the values X a 
which is the value of X in configuration a: 

{X)=Y,P a X a . (3) 

a 

From this follows the rate equation 
d(X) _^dP av 

dt ~ ^ dt a 

a 

= Y, [W a f}P(3 ~ W 0a P a ]X a (4) 

= Y J W a pPp{X a -X l3 ). 

a/3 



In the initial state of the system all the sites are va- 
cant (no particles in the pipe) as we are interested in the 
behavior of the system towards equilibrium. 

B. Master Equation 

Reaction kinetics is described by a stochastic process. 
Every reaction has a microscopic rate constant associated 
with it that is the probability per unit time that the 
reaction occurs. Stochastic models of physical systems 
can be described by a Master Equation. |3| 



C. Master Equation for infinitely- fast diffusion 

We show that we can simplify the Master Equation re- 
moving the fast reactions. In order to remove fast reac- 
tions we stop distinguishing between configurations that 
can be transformed into each other by the fast reactions. 
We split all configurations into disjoint sets such that 
if Ci is one such a set and a, [3 G C;, then a can be 
transformed into (3, or vice versa by fast reactions. If we 
denote 

TTj = Pa, 

aeCi 




\ adsorption 




FIG. 1: Picture of a Single-File System with two types of adsorbed particles 



we can derive the Master Equation for in 



vector P, we can write the Master Equation as 



d^i \ dP a 
~~dt ' ~ "df 

aeCi 

= E E \ W *pPp - Wp a P a ] 

= E E E ^p p p - w m 
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(5) 



With Wij = J2 a eC z Waff-^. 

We see that all fast reactions have disappeared; they 
only contribute to oju, which can be left out of the ME. In 
order to calculate w, 3 , we need the conditional probabil- 
ities — . Because we have fast reactions connecting the 
/3's in Cj we may assume that these d's are in steady-state 
with respect to each other. Hence, the conditional prob- 
ability ^ is nothing but the probability of 3 in steady- 
state if we restrict ourselves to the configurations in Cj. 



D. Simulation methods 



where R and W are time independent. We also introduce 
a new matrix Q, Q(t) = exp[— Ri]. 

This matrix is time dependent by definition and we 
can rewrite the Master Equation in the integral form 



P(i) = Q0)P(0) + / dt'Q{t - f')WP(t'). (7) 
Jo 

By substitution we get from the right-hand-side for P(t') 



P(t) = [Q(t) 

+ f dt'Q(t-t')WQ(t') 
Jo 

+ [ dt' [ dt"Q(t - t')WQ(t' - i")WQ(i") 

Jo Jo 

+ ...]P(0). 



(8) 



1. Dynamic Monte Carlo 

DMC methods allow us to simulate the system gov- 
erned by the Master Equation over time. We simplify the 
notation of the Master Equation by defining a matrix W 
containing the rate constants W a », and a diagonal ma- 
trix R by R a p = ^ 7 W-yp, if a = ft, and otherwise. 
If we put the probabilities of the configurations P a in a the probability distribution 1 — exp[—R aa i\. At time t' a 



Suppose at t = the system is in configuration a 
with probability P a (0). The probability that, at time 
t, the system is still in configuration a is given by 
Qaa(t)P a {0) = exp(— R aa t)P a (0). This shows that the 
first term represents the contribution to the probabilities 
when no reaction takes place up to time t. The matrix W 
determines how the probabilities change when a reaction 
takes place. The second term represents the contribution 
to the probabilities when no reaction takes place between 
times and t', some reaction takes place at time t', and 
then no reaction takes place between t' and t. The sub- 
sequent terms represent contributions when two, three, 
four, etc. reactions take place. The idea of the DMC 
method is not to compute probabilities P a (t) explicitly, 
but to start with some particular configuration, represen- 
tative for the initial state of the experiment one wants to 
simulate, and then generate a sequence of other configu- 
rations with the correct probability. The method gener- 
ates a time t' when the first reaction occurs according to 



reaction takes place such that a new configuration a' is 
generated by picking it out of all possible new configura- 
tions (3 with a probability proportional to W a > a - At this 
point we can proceed by repeating the previous steps, 
drawing again a time for a new reaction and a new con- 
figuration. [Q, [j| One of the most popular DMC method 
in the literature is Random Selection Method (RSM) j|. 
We use this method to simulate the Master Equation of 
our system. 



2. A Dynamic Monte Carlo algorithm for infinitely-fast 
diffusion 

In section C we have derived the ME distinguishing 
between configurations that can be transformed into each 
other by fast reactions. Starting from the ME(5) we give 
a DMC algorithm for simulating infinitely-fast diffusion. 
For our model, diffusion is much faster than adsorption 
and desorption, so the sets are all configurations that 
are connected by diffusion. These sets can be labeled 
only by the sequences of particles A and B, as only the 
number of A and B particles and their order in the pipe 
is important to distinguish the configurations within a 
set. Moreover, all probabilities of configurations within 
a set are the same. This means that ^f- is the number of 

configurations in Cj . There are J ways to distribute 
n particles over S sites. We have then 



Tt-n 



with Pfj £ C n . The summation X^/sec sums over 

configurations. However, for adsorption the left-most 
or the right-most site should be vacant. This gives us 
S — 1 \ 

possible configurations. Each of these gives just 



one configuration in the summation over a. So we get 



S-l 

n ) S - v. 

toads = 2W ad s v ,„ N 7 = 2W ads ^ r - = 2W ads (l - 6), 



S 



.(9) 

where u> ads is the transition probability for the transition 
from a configuration within the set with n particles to a 
configuration within the set with n + 1 particles. Similary 
we find that 



UJdes = 2Wde 



S-l\ 

n-l) 2W des n 



2W des 0. (10) 



The Dynamic Monte Carlo (DMC) algorithm that we 
have used to simulate the system consists of the following 
steps: 

1. Compute the time for the next adsorption or des- 
orption. If the current time is t then the time for that 
process is t + At with 



At 



In i 



2(1 - 9)W ads + 26W dc 



(11) 



where r is a random number picked from the uniform 
distribution on the interval (0, 1] and 9 is the probability 
that the marginal site is occupied. With infinitely fast 
diffusion this probability is given by 9 — (Na + Nb)/S 
with Na and Nb the number of A's respectively B's in 
the system. 

2. Compute for each A in the system a time when it will 
transform into a B. This time is given by t + r with 



In 7' 



(12) 



In this expression P stands for the probability that the 
A is at a protonic site. If we number the particles in the 
system from left to right 1,2,..., Na + Nb and the sites 
also from left to right 1, 2, . . . , S then the probability that 
particle number n is at site number s, P*, is given by 



P.; 



s-lW S-s 
n-l) \N A + N b - n 

S 

Na + Nb 



(13) 



P for particle n is then the sum of this expression over 
all protonic sites 



s=l 



(14) 



In both expressions 9 = ft is the coverage. 



where 8 S — 1 if site s is protonic. 

3. Change those A's with r < At into B's. 

4. Determine the next process at the marginal sites. 
It is an adsorption with probability proportional to 
(1 — 9)W a ,ds and a desorption with probability propor- 
tional to 6Wdcs- The process is equally likely to occur on 
the left- or the right-hand-side. 

5. Change the number of particles in the system accord- 
ing to the next process at the marginal sites. 

6. Update the time. 

7. Repeat steps 1 to 6. 



III. ANALYTICAL RESULTS 

In this section, for some special cases such as low load- 
ing limit and fast and slow reaction, we are able to derive 
some expressions for the productivity in the steady-state. 
For the case of low loading limit we can also derive the 
rate equations of the system. 



A. The low loading limit. 

In this case we can assume that there is never more 
than one particle in the system. The following rate equa- 
tions then hold. 

d(A) _ 2 Wa ds _ ^ _ 2W dcs ^ 



dt 



S 
W r> 



T <A> 



(15) 



W, 



d(B) = 2W dcs 

dt S ' 27Viner, 



T (A> 



Here (X) is the probability that there is an X at an ar- 
bitrary site. For steady state we get 



(A) 



2W ads W dc s 



(Wads + W d cs)(2W d cs + JVprotWrx) 



NpiatWedsWrx 



(16) 



(Wads + W dc s)(2Wdes + ^pro t W rx ) 

From this we immediately get the turnover frequency 

2WadsWdesW rx 



WTOF 



(W ads + W d es)(2W dcs + iV prot W rx ) 



(17) 



We see that in this limit the turnover frequency does not 
depend on the system size. 

Comparing the number of B's produced from the an- 
alytical results with the DMC results for the case when 
W ads = 0.0033, Wdos = 0.9967, W rx = 0.1 and different 
distributions of the reactive sites, we remark that we get 
similar results. 



B. Fast and slow A — > B reaction. 

If the reaction is fast, and there are not too many par- 
ticles in the system, then all particles in the system are 
Bs. This means 



W 



ads 



and 



WTOF 



Wads + Wies ' 



1 2WadsW dc s 
Nprot Wads + W d e 



(18) 



(19) 



The restriction of not too many particles is necessary, 
because particles should all always be able to reach a 
protonic site. This means 



W 



ads 



W 



ads 



w dc 



S 



(20) 



must hold. This relation depends on the distribution of 
the reactive sites. When reaction is fast this means that 
it depends on the distance from the margins to the first 
protonic site. 



Comparing the site occupancy with B's(20) with re- 
sults from the DMC simulations, we find that for the 
reaction rate constants (W a ds=0. 03333, Wd os =0. 96667, 
W rx =10 and all the sites reactive, the results are similar. 

If the reaction is slow, then there are only occasionally 
B's in the system. This means 



(A) 



W 



ads 



Wads + W dc 



(21) 



All particles in the system will be renewed between two 
subsequent formations of a B. Therefore 



WTOF 



W a dsW rx 
Wads + W dc 



(22) 



Comparing the site occupancy with A's(23) with re- 
sults from the DMC simulations, we find that for the 
reaction rates constants W a ds=0. 03333, Wd C s=0. 96667, 
W rx =0.001 and all the sites reactive, the results are sim- 
ilar. 



IV. SIMULATION RESULTS AND 
DISCUSSION 

A. All sites reactive 

We remark that DMC methods with regular high rates 
for diffusion tend to give similar results as DMC method 
for infinitely- fast diffusion described in section |IID 1 . 



The results of this comparisons are in figure |[ We con- 
clude that the DMC method for infinitely-fast diffusion 
we have defined is a correct method to simulate the be- 
havior of the system in the limit W d \s —> oo. 

In Q we have seen that for the case when all the sites 
are reactive, the site occupancies of the system obtained 
from DMC simulations show that the system is not ho- 
mogeneous even for very fast diffusion rates. We find 
the same effect also using DMC for infinitely-fast diffu- 
sion, for different loadings and for different reaction rate 
constants. 

We study also the dependence of the occupancy profiles 
on the reaction rate constant at different loadings of the 
system. 

The simulation results in figure ^| show how the shape 
of the profiles changes with reaction rate W rx for high 
and low loading of the pipe, when all the sites are reac- 
tive. We find that for high loadings, as an effect of the 
blocking, the middle sites have the same probability to 
be occupied for both fast and slow reactive system. Only 
the occupancy of marginal sites is influenced by the reac- 
tivity, such that in fast reactive systems we have a higher 
probability to have a B near the marginal sites, and, in 
consequence, a better productivity. For slow reactive sys- 
tems, the occupancy profiles are scaled with reaction rate 
constant. We notice that the productivity is growing in- 
creasing reaction rate in case of low loading systems al- 
most as fast as in the case of the high loading systems 
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FIG. 2: a)Dynamic Monte Carlo results for site occupancy 
for the case of infinitely-fast diffusion and low loading ( W a ds = 
0.2, Wdes ~ 0.8) of a system of length S = 30. The continuous 
line is for the site occupancy for W rx = 0.1, the first dotted 
line near the continous line is for W rx = 0.2, and the second 
is for W rx = 0.4. b)Dynamic Monte Carlo results for site 
occupancy for the case of infinitely-fast diffusion and high 
loading (Wads = 0.8, Wdes = 0.2) of a system of length S = 
30. The continuous line is for the site occupancy for W rx = 
0.1, the first dotted line near the continous line is for W rx = 
0.2, and the second is for W rx = 0.4. 



FIG. 3: Profile occupancies for a system of length S — 30, 
Wads=0.6, Wdes=0.4 and W rx =0.1. The continous lines are 
the profile occupancies for A (the lower) and B (the higher) 
using DMC for infinitely-fast diffusion. The dotted lines are 
the profile occupancies for A(the lower) and B (the higher) 
using DMC with a regular high rate for diffusion (Wdifr = 
1600). 
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because of the diminished effect of the blocking in the 
pipe. Comparing for instance the rate o growth for B 
production when reaction rate constant is increasing from 
W rx = 0.1 to W rx = 0.4 " B%f° A )/S), in 

case of low loading (6 = 0.2) and high loading (6 = 0.8), 
we find almost the same rate of growth in both the cases, 
and this is 0.5. 

In figure || we have the logarithmic shape of the oc- 
cupancy profiles for A and B. These profiles show that 
smaller systems are less reactive as less As are inside the 
pipe. This explains the faster decrease of (A„) for small 
systems in figure |J. For infinitely-long systems we expect 
to have straight lines corresponding to an exponential de- 
crease of (An). 



B. Only some of the sites reactive 

For all the sites reactive, in paper [El, we have shown 
that even when Waff — > oo, DMC results with regu- 
lar high diffusion rates indicate that the system doesn't 
become homogeneous. Using DMC for infinitely-fast dif- 
fusion we find that also for different distribution of the 



FIG. 4: The logarithm of the DMC(Random Selection 
Method) profile occupancies ((A„)) for W a ds = 0.6, Wdes = 
0.4, Wrx = 0.1, and when all the sites are reactive, for various 
system sizes S. 



reactive sites we find a non-homogeneous distribution of 
the particles in the system. In figure || we can see that 
for marginal sites reactive and for middle sites reactive 
as well, we have single-file effects also for infinitely-fast 
diffusion for different rates of reactions. 

We notice that the profiles in case of high loadings are 
very slowly dependent on conversion, both for middle and 
for marginal sites reactive. Comparing with simulation 
results for occupancy profiles in jj], where middle sites 
and marginal sites are reactive, and for regular diffusion 
rate constants in the domain (2 . . . 10), we notice that the 
profiles are similar. We can conclude that as the effect of 
infinitely-fast diffusion is absent for Single-File Systems, 
the diffusion is not so important in these systems. 

For low loading, when marginal sites are reactive, the 
occupancy profiles are scaled with W rx - We notice that 
conversion in figure]^ is the rate determining step. In this 
case, the middle sites doesn't have the same occupancy 



for different reaction rate constants like in the case of all 
the sites reactive but are strongly dependent on W rx . 

We notice also that the productivity in case marginal 
sites are reactive is growing faster increasing the reac- 
tion rate constants at low loadings than at high loadings. 
Comparing, for instance, the rate o growth for B pro- 
duction Bp ro d when reaction rate constant is increasing 
from W rx = 0.1 to W rx = 0.4, in case of low loading 
(0 = 0.2) and high loading (6 = 0.8), wc find that the 
rate of growth of B productivity at low loadings (0.3) is 
larger than at high loadings (0.25). 

When we have middle sites reactive we have higher 
probability to find A's near marginal sites. The produc- 
tivity is smaller than in all the other cases. The profiles 
are again scaled with W IK for low loadings. For high 
loadings and middle sites reactive, the profiles for differ- 
ent conversion rates are almost similar, so the productiv- 
ity can only be increased in this case only increasing the 
number of reactive sites. 



V. SUMMARY 

We have used DMC methods and analytical tech- 
niques to analyze Single-File Systems for which diffusion 



is infinitely- fast. We simplified the ME removing fast 
diffusion and we have presented a DMC algorithm for 
infinitely-fast diffusion that simulate this ME. We show 
that DMC with regular high rates gives the same results 
as DMC for infinitely-fast diffusion. The fundamental 
assumption considered for infinitely-fast diffusion in the 
analytical results is that all configurations related by dif- 
fusion have the same probability. 

In the limiting cases such as low loading limit and slow 
and fast conversion, we are able to derive expressions for 
the B productivity. We notice that the number of Bs 
produced per unit time in these cases doesn't depend 
on the system size. Comparisons between analytical and 
DMC results reveal similar results for the productivity. 

DMC results show that when all the sites are reactive 
and when only some of the sites are reactive, diffusion 
has no influence on the single-file properties of the sys- 
tem. Different results for the dependencies of the oc- 
cupancy profiles and productivity on the reaction rate 
constant and different distributions are categorized. The 
occupancy profiles show that smaller systems are less re- 
active as less As are inside the pipe. 
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